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INTRODUCTION 


^ Under  NRL  contract  #N00173-77-C-0215,  JAYCOR  has  performed  an  analytic  and 
experimental  study  of  some  of  the  fundamental  aspects  of  cyclotron  masers 
(gyrotrons).  In  particular,  JAYCOR  has  designed  and  tested  a high  power  gyrotron 
oscillator,  and  has  designed  and  used  a gyrotron  testing  facility.  This  report 
discusses  the  results  of  this  project  in  detail. 

The  JAYCOR-NRL  contract  includes  the  following  work  statement: 

(1)  Design  and  test  a 200  kw  cyclotron  maser  testing  facility  which 
includes  the  modulator,  pulse  transformer,  power  supply,  superconducting  magnet, 
and  auxiliary  subsystems. 

(2)  Design  and  test  maser  components  to  include  the  RF  interaction  structure 
and  output  coupling  network. 

(3)  Design  and  test  the  RF  instrumentation  package  to  include  a millimeter 

wave  calorimetric  load  and  an  output  mode  analyzer, 

(4)  Teat  and  evaluate  the  magnetron  injection  gun  to  be  used  in  the 
cyclotron  maser.  Compare  measured  performance  with  theoretical  predictions 

and  optimize  gun  performance. 

(5)  Design  and  test  the  RF  driver  package. 

Essentially  this  work  statement  calls  for  the  design  of  a laboratory  for 
the  testing  of  cyclotron  masers  and  the  design  and  test  of  the  various  components 
of  an  initial  device.  One  element  not  defined  in  the  above  statement  was 
whether  the  initial  4syice  was  to  be  an  oscillator  (gyronmonotron)  or  amplifier 
(gyro-TWA) . After  some  initial  study  and  consultation  with  NRL  scientists,  it 
was  decided  that  an  oscillator  was  the  best  Initial  device,  in  that  it  was  the 
simpler  of  the  t a,  and  most  likely  to  quickly  yield  the  most  information  about 
the  basic  physical  mechanisms.  No  RF  driver  was  needed  for  an  oscillator, 
therefore  this  report  will  not  deal  with  one. 


While  preliminary  plans  for  a system  to  analyze  the  electron  beam 
characteristics  were  drawn  up  under  this  contract,  it  was  determined  by  NRL 
scientists  that  it  was  not  feasible  to  construct  the  device  at  that  time. 
Therefore,  no  analysis  on  the  beam  characteristics  was  carried  out.  Instead, 
more  extensive  testing  of  the  gyrotron  system  as  a whole  was  done  than  was 
called  for  in  the  contract.  This  work  is  reported  in  Section  III. 

The  body  of  this  report  will  be  broken  down  into  two  sections,  and  an 
appendix.  The  first  section  is  a short  discussion  of  the  linear  theory,  the 
second  a report  on  the  nonlinear  analysis  and  design  work.  The  appendix 
describes  the  gyrotron  test  facility  which  has  been  established  at  NRL. 
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I.  Basic  Mechanism  and  Linear  Calculation 
Basic  Mechanism 

The  gyrotron,  or  cyclotron  resonance  maser,  mechanism,  is  a relativistic 
Instability  involving  a fast  waveguide  mode  and  the  fast  cyclotron  mode. 
These  modes  are  given  by 


<jj2  - k2c2+  wco2 

1.1 

u - kv(|  + Slc/ y 

1.2 

respectively,  w is  the  wave  frequency,  in  radians  per  second,  k is  the  wave- 
number,  c the  speed  of  light,  ojco  the  cutoff  frequency  of  the  waveguide, 
vw  the  parallel  velocity  of  the  beam,  y the  relativistic  mass  factor,  and  &c 
the  classical  cyclotron  frequency.  These  dispersion  relations  are  shown 
graphically  in  Figure  1.1a.  The  modes  are  shown  uncoupled. 

The  actual  mechanism  is  discussed  extensively  by  Sprangle^^.  In  summary, 
the  instability  results  from  electron  bunching  that  occurs  because  the  relati- 
vistic cyclotron  frequency  is  inversely  proportional  to  y.  The  energy  to 
drive  the  instability  comes  from  the  motion  of  the  electrons  perpendicular  to 
the  magnetic  field,  which  couples  to  the  wave  electric  field  in  a manner  shown 
in  Figure  1.1b.  The  electrons  in  the  figure  are  shown  already  bunched — the 
bunching  mechanism  is  discussed  by  Sprangle. 

The  geometry  of  the  model  used  for  both  the  linear  and  nonlinear  calcula- 
tions is  shown  in  Figure  1.2. 

Linear  Calculations 

The  linear  theory  of  gyrotron  oscillators  with  right  circular  cavities  has 
been  extensively  developed  by  Chu(2)  and  Sprangle^.  Chu's  work  will  be 
extensively  referenced  in  this  report.  For  the  sake  of  completeness,  a summary 
of  their  results  which  are  used  in  the  design  methods  of  this  report  will  be 
given. 
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Dispersion  curves  of  the  waveguide  (upper)  and  fast 
cyclotron  modes  (lower) 


Mechanism  of  energy  transfer  from  the  electrons  (dots)  of  a 
bunched  annular  electron  beam  to  the  electromagnetic  field  via 
the  azimuthal  electric  field.  This  particular  diagram  is  for 
the  TE011;  u)£sftc.  Note  that  for  this  interaction  the  energy 
transfer  is  essentially  continuous,  as  the  phase  of  the  Eg  field 
shifts  with  the  change  in  electron  velocity. 


Figure  1.1b 


Li  * 

.ci<r  ^ 

An  essential  point  in  the  linear  theory  is  the  coupling  coefficient,  F, 
which  gives  the  rate  at  which  energy  is  transferred  from  the  wave  to  the 
electrons.  The  coefficient  is  a function  of  the  cavity  and  beam  geometry, 
the  beam  energy,  current  and  velocity  distribution,  the  waveguide  mode,  the 
harmonic  number  of  the  cyclotron  frequency,  and  the  phase  shift  between  the 
beam  and  waveguide  modes.  For  clarity,  we  define  these  parameters  below,  and 
discuss  them  briefly. 

b ■ cavity  radius 

L • cavity  length 

r ■ average  beam  radius 
b 

V • beam  voltage 

6 ■ beam  velocity  normalized  to  the  speed  of  light  in  vacuum,  c 


Pjj  * beam  power 
Y ■ relativistic  factor  ■ 


m,n,£  ■ radial,  azimuthal,  and  axial  eigen-numbers  of  a cylindrical  cavity 
s * harmonic  number  of  the  magnetic  field 

A ■ (to  - kv  - 1.3 

" y 

The  parameter  A is  the  phase  shift  between  the  waveguide  and  beam  modes,  u>  is 

given  by  (1.1),  k the  wave  number,  is  given  by 

, JLtt 
L 

The  parametrical  variations  in  Chu's  work  that  are  important  in  the 

design  methods  discussed  here  are  those  for  the  coupling  coefficient  as  a 

function  of  the  beam  radius,  and  the  phase  shift.  The  coupling  coefficient 

is  a maximum  for  certain  beam  positions;  the  beam  positions  required  for  the 

(2 ) 

maxima  are  given  in  the  reference  • For  the  TEg^mode,  with  s«l, 
the  optimum  beam  position  is  r^  ■ 0.47  b. 
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The  dependence  on  the  phase  shift  is  most  usefully  illustrated  through 
a plot  of  the  threshold  condition.  The  threshold  condition  is  defined  in 
terms  of  the  minimum  beam  power  and  cavity  Q required  for  oscillation,  i.e., 

V>  i <pbQ>th  !•* 

and  is  inversely  proportional  to  the  coupling  coefficient.  A plot  of  the 
threshold  (P^Q)  for  a cavity  used  in  the  35  GHz  oscillator  is  shown  in  Figure 
1.3. 


II.  Non-Linear  Calculations 

While  an  estimate  of  the  maximum  achievable  efficiency  of  an  oscillator 
is  possible  analytically,  an  accurate  quantitative  determination  of  this 
parameter  has  had  to  be  done  numerically.  A code  for  determining  the  energy 

loss  of  charged  particles  due  to  the  interaction  with  guide  electromagnetic 

(1) 

fields  has  been  developed  by  Drobot  » and  modified  for  use  in  cavity 
calculations. 

This  code  utilizes  a routine  "CYCMASC"  which  integrates  the  electron 
energy  loss  due  to  the  interaction  with  the  RF  fields  along  a self  consistent 
orbit.  The  code  is  2-dimensional,  modeling  only  azimuthally  symmetrical 
fields  and  particle  distributions. 

Fully  relativistic  equations  were  used  in  the  code.  For  these  calcula- 
tions the  average  amplitude  of  the  electromagnetic  field  is  held  constant 
(in  time);  i.e.,  the  system  is  assumed  to  be  in  a steady  state.  The  steady 
state  condition  is  made  self  consistent  by  assuming  that  all  of  the  power  lost 
by  the  electrons  is  transferred  to  electromagnetic  fields  and  radiated,  i.e., 
the  system  is  lossless.  The  calculation  can  readily  be  altered  to  include 
cavity  losses,  which  must  be  considered  for  very  high  frequency  or  low  power 


devices. 


Plot  of  the  threshold  value  for  QP^.  in  kw,  versus  the  phase 
shift  A,  for  the  TEq..  mode  in  a cavity  with  L/b  = 5.4,  and  a 
beam  with  Y ■ 1.14  for  •<*  1.5. 


Inputs  to  the  Code 

A simplified  flow  diagram  of  the  code  is  given  in  Figure  2.1. 

The  inputs  to  the  program  are:  the  initial  electron  parameters  a and  y, 
the  electron  beam  current,  the  cavity  mode  numbers  n and  £ (m-0  always)  the 
shape  and  magnitude  of  the  RF  electric  field,  the  magnitude  and  shape  of  the 
DC  magnetic  field,  and  the  cyclotron  harmonic  number  s.  We  will  discuss 
these  parameters  in  detail. 

The  ratio  of  perpendicular  to  parallel  electron  momenta  a,  is,  as  in 
the  linear  theory,  chosen  from  electron  gun  requirements.  Clearly  it  is 
desirable  for  a to  be  as  large  as  possible.  However,  present  gun  designs 
seem  limited  to  maximum  a's  of  1.5  to  2.0.  All  of  the  code  runs  reported  here 
are  for  a equal  to  1.5.  The  choice  of  the  total  electron  energy,  given  by 
Y , is  not  so  limited.  A lower  limit  to  Yq»  may  be  arrived  at  through  the 
consideration  of  the  required  product  of  beam  power  and  cavity  Q for  the  initia- 
tion of  oscillation,  derived  from  the  linear  theory.  This  product  decreases 
slowly  with  decreasing  beam  voltage,  so  that,  with  constant  Q,  increasing  beam 
currents  are  required  with  decreasing  Yq.  Eventually  the  required  current 
will  exceed  that  which  can  be  propagated,  and  prior  to  that  significant  space 
charge  effects  may  be  expected.  Chu^)  has  given  an  estimate  for  an  upper 
bound  for  a space  charge  when  it  is  manifested  as  spreads  in  the  steady  state 
electron  energy  and  parallel  velocities.  This  limit,  in  terms  of  a spread  in 
the  phase  shift  A,  is 


6A  ■ (£7r68  + s 


5* 

c 


1 £u 

a TT7 


1 

* 


«2.5tt 
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The  magnitude  of  the  electromagnetic  field  necessary  may  be  estimated 


from  the  requirements  that 


dl  - (Y„-l)  — 

o e 


where  the  integration  is  along  an  electron  orbit  and  the  effective 

electric  field,  is  given  by  (2) 


/l6Pt V 

I FLb2  Jq2  (x^) J 


where  L is  the  cavity  length,  Jq  (xn)  is  the  Bessel  function  of  the  first 
kind,  and  T is  the  electron  transit  time.  The  coupling  coefficient  F is 
given  by  Chu(2)t 

The  magnitude  of  the  magnetic  field  can  again  be  estimated  from  the  linear 
theory.  For  weak  interactions  (*J«0.1)  the  most  efficient  operation  will  be 
obtained  where  the  beam-wave  coupling  is  maximum.  This  corersponds  to  the 
function  F being  the  most  negative,  by  definition.  This  will  occur  in  general 
for  the  phase  shift  A>0  where  A is  defined  by  equation  (1.3).  For  strong 
interactions,  (y  >.l)  the  electrons  will  lose  energy  as  they  traverse  the 
cavity.  For  the  parameter  range  of  interest  A will  decrease  with  Y«  Assuming 
that  the  axial  magnetic  field  is  constant,  an  estimate  for  the  initial  A that 
is  slightly  higher  than  that  for  optimum  coupling,  thereby  allowing  for  some 
decrease  in  A,  seems  reasonable.  A more  effective  way  to  optimize  the  inter- 
action would  be  to  axially  vary  the  magnetic  field,  thus  holding  A constant. 
This  will  be  the  subject  of  a future  investigation. 

The  choice  of  the  harmonic  number  s i3  dependent  upon  the  available 
magnetic  field.  In  most  applications  it  is  expected  that  it  will  be  desirable 
to  use  as  low  a field  as  possible.  This  must  be  balanced  against  the  diffi- 
culties caused  by  the  Increasing  density  of  competing  modes. 


e 
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Outputs  from  the  Code 

The  outputs  of  the  code  are  the  efficiency  of  energy  extraction  n,  and 
the  required  Q of  the  cavity.  The  efficiency  is  defined  as 


- (XL) 


where  y^  and  y^  are  the  initial  and  final  relativistic  velocity  factors, 
respectively. 


The  cavity  Q required  to  achieve  the  given  efficiency  is  calculated  (in 
the  code)  using  the  definition 


where  W^  is  the  stored  electromagnetic  field  energy.  For  a right  circular 


cavity  and  a TE  mode, 

oni. 

v - r 


[J‘  T1 '] 


where  p is  the  n n zero  of  the  bessel  function  J and  Z is  the  characteristic 
on 

impedance  of  the  cavity. 

We  note  that  there  is  no  guaranty  that  the  required  Q can  be  realized 

in  practice.  In  fact  the  realizable  Q has  both  an  upper  and  lower  bound. 

( 3) 

The  minimum  Q achievable  is  usually  termed  the  diffraction  limited  Q , 
and  is  given  by 


■JV- 


« A2L  (A  )2  (a >«<*>  ) 
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This  is  the  Q of  a smooth  waveguide  of  length  L with  no  reflection  at  either 
end.  For  most  cavities  considered  for  gyromonotrons,  L/^>5.  For  a one-half 
guide  wavelength  cavity  (TEon^)  the  diffraction  Q is  then  >,300.  We  note 
that  since  the  area  integrals  in  (2.9)  cancel,  the  diffraction  Q is  independent 
of  the  radial  eigen-number  (n). 

An  upper  bound  for  the  cavity  Q is  given  by  the  losses  in  the  cavity 

walls,  i.e.,  by  the  ohmic  Q.  For  a right  circular  cavity  with  a TE  mode 

onJt 


Q 


ohmic 


wpa  jkQn  + 

2R  (k*n  + 2 ttH2  <£)*J 


2.10 


where  R is  the  surface  resistance  of  the  cavity  walls. 

With  copper  walls  the  expression  becomes 

15.2  ff  b (k?on  + TT2*2  (b)2 

Wc  * “ L_  2.11 

<kon  + 2lr2*2  (£)3 

where  kon  is  the  nch  non-vanishing  2ero  of  (kon) » Bessel  function  of 
the  first  kind.  At  35  GHz  a TEqii  cavity  with  a radius  of  0.5  cm  and  (b/L) 
=0.2  will  have  an  ohmic  Q of  14,000. 


Results  of  the  Simulations 

It  is  clear  from  the  formulation  of  the  cavity  problem  that  the  gyro- 
monotron  efficiency  is  a function  of  many  parameters.  Specifically,  the 
efficiency  is  dependent  upon  b/L,  A,  E^,  rb/{j  » Yo*  T»  an<*  Eo>  an<*  Yo*  As 
We  have  examined  numerically  the  effect  of  variations  in  b/L,  E^,  andY*  As 
mentioned  above,  rb/.  la  chosen  analytically.  °<  is  assumed  to  be  1.5,  a 

D 

value  deemed  reasonable  from  present  electron  gun  designs. 


.1 
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Usually  32  particles  were  used  in  the  simulation.  Larger  numbers  of 
particles  were  occasionally  run  as  a check,  but  no  significant  change  in  the 
results  was  observed.  The  particles  were  initially  uniformly  distributed  in 
phase,  but  injected  into  the  cavity  at  the  same  time.  Again,  a check  was 
occasionally  made  on  this  simplified  condition  by  Injecting  several  groups 
of  32  particles  at  different  times.  No  significant  change  in  the  code  output 
was  observed.  The  modeling  of  the  beam  electrons  being  initially  randomly 
distributed  in  phase  is  justified  since  the  beam  is  well  randbmized  by  the 


time  it  reaches  the  cavity. 


(5) 


The  code  was  checked  by  comparing  its  results  with  those  of  the  linear 
theory  in  the  regime  of  low  Eg  (low  QP),  where  the  linear  theory  is  valid. 
Agreement  was  excellent,  as  can  be  seen  in  Figure  2.2. 

Output  from  the  code  was  (in  part)  in  the  form  of  Figure  2.3.  This 
particular  output  is  for  a case  where  strong  phase  trapping  was  observed. 

The  large  curve  is  a plot  of  efficiency  vs.  axial  position.  The  inserts  are 
plots  of  the  particle  perpendicular  noraallized  velocity  (8^)  versus  phase. 

As  mentioned  above,  the  particles  are  injected  with  distributed  phase  and 
uniform  8j_.  Interaction  with  the  cavity  field  causes  bunching  and  energy 
extraction  from  the  electrons.  (Figures  2.3,  b,  c).  The  maximum  energy 
extraction  occurs  at  the  axial  position  m . At  this  point  the 

electrons  are  well  phased  trapped.  Some  particles  regain  energy  as  they  drift 
to  the  end  of  the  cavity  (Figures  2.3,d). 


t: 


Variations  of  Parameters 

The  product  of  the  cavity  Q and  the  beam  power  is  proportional  to 
a) J* Eg2  2TrrdrdA  and  is  therefore,  for  a given  cavity  geometry  and  efficiency. 


a measure  of  E, 


A plot  of  r|  versus  QP  for  b ■ 0.53  cm  (cavity  resonant 
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Figure  2.3.  Computer  output,  showing,  in  the  large  plot,  the  gyrotron  efficiency 
versus  the  number  of  cyclotron  periods  an  electron  has  undergone 
after  entering  the  cavity.  The  end  of  the  cavity  corresponds  to  the 
right  hand  limit  of  the  plot.  Note  that  the  efficiency  reaches  a 
maximum  at  Z « .75  L,  then  drops  due  to  phase  trapping.  The  small 
plots  show  the  perpendicular  velocity  (versus  azimuthal  phase)  for 
three  different  times.  At  t ~ 0,  the  electrons  are  unbunched,  and 
all  have  the  same  perpendicular  energy.  At  t * 2.13,  the  electrons 
are  bunched,  and  most  have  lost  energy,  and  at  t * 2.94,  the  bunch 
has  started  to  trap. 


frequency  equal  to  35  GHz)  for  a TEqu  mode  is  given  in  Figure  2.4.  These 

curves  are  for  cases  where  the  parameter  A has  been  adjusted  for  optimum 

efficiency.  It  is  observed  that  an  optimum  occurs  for  each  L/b,  with  these 

optima  increasing  with  this  ratio.  The  Q required  for  optimum  efficiency 

decreases  EgdU  with  L/b.  This  is  consistent  with  the  concept  that  there  is 

an  optimum  E.df  , below  which  not  all  available  energy  has  been  extracted 
0 

from  the  electrons,  and  above  which  phase  trapping  occurs.  The  effect  of 

L L 

trapping  can  be  seen  in  the  cyclic  nature  of  the  curves  for  — - 11  and  ~ * 7.8. 
A corresponding  bounce  is  observed  on  the  n vs.  z curve  for  simulation  runs 
for  high  values  of  QP. 

n vs  A 

It  is  expected  that  for  weak  (i.e.,  low  efficiency)  interactions  the 
optimum  initial  A should  be  equal  to  that  giving  the  strongest  linear  coupling 
coefficient.  In  these  interactions  A will  remain  essentially  constant  as  the 
electron  traverses  the  cavity,  since  y remains  almost  unchanged  (see  equation 
1.3).  For  higher  efficiency  interactions  there  will  be  a significant  drop 
in  y s an  electron  crosses  the  cavity,  and  a corresponding  decrease  in  A. 

In  order  to  maintain  A in  the  range  where  strong  negative  coupling  results 
for  the  maximum  time  it  is  necessary  to  increase  the  initial  phase  shift 
from  the  optimum  linearly  calculated  value  (see  Figure  2.5). 

n vs  y 

An  estimate  for  the  dependence  of  the  efficiency  on  the  beam  energy  can 


be  derived  by  noting  that  energy  can  be  extracted  from  the  electrons  while 

(2) 

-0.5Tr<A<2ir.  The  resulting  expression  for  the  efficiency, 


EFFICIENCY 


Figure  2.5.  Plot  of  the  electronic  efficiency  versus  the  phase  shift  A 
and  the  corresponding  detuning  in  magnetic  field,  given  as 
over  the  operating  frequency  u>.  TEqj^,  L/b  * 5.0, 

QP.  - 250,  Y - 1.14,  <*-  1.5. 


2- 5Try0<Yf  > 
s ft  X (Y„-l) 


is  only  an  upper  bound,  and  depends  on  one's  estimate  for  the  average  of 
electron  energy  at  the  end  of  the  cavity  (<y^>).  In  many  designs  y o and 
<Y^>-1.  Here  the  dominant  variation  in  n will  come  from  the  term  (Yo-1) . 

In  this  case  n decreases  with  increasing  Yo • Simulation  results  qualitatively 
confirm  this  dependence,  as  is  shown  in  Figure  2.6. 


H vs.  s 

Equation  2.9  predicts  a constant  efficiency  for  constant  s ftc.  There- 
fore, ignoring  the  effects  of  temperature,  spaee  charge,  and  competing  modes, 
it  would  appear  that  operating  at  harmonics  of  the  cyclotron  frequency  causes 
no  reduction  in  efficiency.  This  appears  to  be  essentially  correct  at  higher 
voltages.  (>_  30  kv)  as  is  shown  in  figure  2.6.  However,  it  is  expected  that 
space  charge  and  temperature  effects  will  become  more  important  with  increas- 
ing s. 

Field  Shape 

The  shape  of  Eq,  as  a function  of  z,  can  be  entered  into  the  code,  either 
point  by  point,  or  by  a function.  As  mentioned  above.  Eg  was  normally  taken  as 

E - E sin  . 

e e0  L 

However,  for  one  cavity  (L/b  ■ 7.8),  the  field  measured  for  a slightly 
tapered  version  was  used  for  E.,  and  the  efficiency  found.  The  field  profile 
is  shown  in  Figure  2.16.  The  maximum  efficiency  found  was  37%,  versus  34Z 
for  the  equivalent  untapered  cavity.  Since  no  other  profiles  were  run,  it 
is  not  known  whether  the  electric  field  profile  chosen  was  an  optimum  one, 
and  perhaps  even  larger  increases  in  efficiency  might  be  obtained  with 
different  profiles. 
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Figure  2.6.  Plot  of  the  efficiency  versus  beam  energy,  for  two  values  of  the 
cyclotron  harmonic  numbers.  Mode  * TEq^,  L/b  * 8,  1.15; 

QP  and  A are  optimized  for  each  point. 


The  reason  for  the  increase  in  efficiency  is  not  all  together  clear. 

A possible  explanation  is  that  the  slow  rise  of  the  field  at  the  entrance 
of  the  cavity  promotes  better  (more  adiabatic)  bunching,  and  that  the  more 
abrupt  drop  of  the  field  at  the  output  reduces  phase  trapping, as  can  be 
observed  in  Figure  2.2. 

Design  Procedures 

As  has  been  shown  in  the  above  sections,  the  efficiency  of  a single 
cavity  gyrotron  oscillator  is  a function  of  many  variables.  These  variables 
are  interdependent,  and  are  bounded  by  external  constraints.  We  present  here 
by  illustration  a rough  method  by  which  an  optimum  design  can  be  achieved, 
with  the  various  constraints  discussed  above. 

Figure  2.7  gives  an  outline  of  the  order  of  calculations  necessary  for 
the  design.  The  first  column  contains  the  independent  variables,  of  which 
f,  P,  s,  and  L are  assumed  to  be  determined  by  desired  output  characteris- 
tics (f,  P)  or  given  system  characteristics  (s,  r,  P).  The  parameters  in 
parentheses  are  ones  which  must  be  varied  in  order  to  achieve  an  optimum 
(and  feasible)  design. 

The  second  column  contains  those  quantities  which  can  be  calculated 
by  analytic  methods.  The  third  column  has  conditions  which,  although  not 
final  design  parameters,  must  not  exceed  given  bounds  in  order  for  the  design 
to  be  feasible.  The  fourth  column  contains  those  qualities  which  are  deter- 
mined using  the  code  CYCMASC.  The  last  column  holds  the  final  design  parameters. 

35  GH2  Oscillator  Design 

Parameter  specifications 

The  above  procedures  were  used  to  design  the  35  GHz  oscillator.  For  this 
device,  it  was  intended  that  the  design  incorporate  an  electron  gun  designed 
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Figure  2.7.  Outline  of  the  procedure  used  for  a gyromonotron  design 
text  for  explanation. 


ohmic 


for  a gyrotron  amplifier.  The  beam  parameters  for  this  gun  were 
V - 70  kv 
1 < 10A 
ct  = 1.5 


For  simplicity  of  operation  and  for  ease  of  handling  of  the  output 

radiation,  the  azimuthal  eigen-number  n was  chosen  to  be  1.  The  axial 

eigen-number  was  initially  chosen  to  be  1,  since  it  was  expected  that  the 

TE mode  would  have  a lower  starting  power  than  the  TE„,  . modes. 

011  01(2,3,4...) 

This  point  was  validated  by  the  simulation  runs,  and  is  discussed  below. 

r/b  was  mentioned  in  section  1,  and  is  optimally  0.47.  (For  the  TEqj^ 
case  the  radius  of  0.47b  corresponds  to  the  position  of  the  peak  of  the 
electric  field.) 


Calculation  of  the  space  charge  parameters,  using  equations  2.1,  2.2, 
and  2.3,  yields 


Ay/y  - 2.9K10' 


-3 


A8/8  - 4.13x10 
z z 

sfl 


-2 


(£ir  AS  + — S L6A  ) _JL  - 0.33 


z c y z 0 


zo 


We  note  that  this  last  quantity  is  much  less  than  2.5tt,  thus  meeting  Chu's 
criterion  for  the  effects  of  space  charge  to  be  negligible. 


Using  the  graphs  of  Figures  2.4  and  2.5,  we  find  the  optimal  values  for 

L/b,  A,  and  QP.  It  is  clear  from  Figure  2.4  that  the  efficiency  increases 

monotonically  with  the  ratio  L/b.  Concommitantly , however,  a decreasing  value 

of  QP  is  required  to  achieve  the  optimum  efficiency  for  that  ratio.  For  a 
b 

given  beam  power,  QP^  has  the  lower  bound  of  , where  ia  given 
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by  equation  2.9.  Thus  there  will  be  a finite  optimal  value  for  L/b.  This 
point  is  illustrated  by  Figure  2.8.  For  P equal  to  0.75,  the  optimum  value 
for  L/b  is  5.4,  and  the  optimum  QP^  is  1.5x10®. 

The  maximum  efficiency  yielded  by  the  code  for  these  values  of  input 
parameters  is  27%.  The  magnetic  field  which  yields  this  efficiency  is 
approximately  13.53  kg  and  the  corresponding  phase  shift,  A,  is  2.8. 

The  threshold  value  for  QP,  is  found  to  be  approximately  6x10^.  The 

D 

operating  QP  is  2.5  times  this  value,  and  is  thus  satisfactory  from  this 
standpoint. 

The  final  design  parameters  are  listed  in  table  2.1.  Note  that  the 
gyrotron  will  have  an  output  power  of  150  kw  at  the  peak  efficiency  of  27%. 
The  NRL  electron  gun  is  capable  of  producing  a beam  of  up  to  15A  (for  short 
periods),  thus  an  output  power  of  268  kw  can  be  achieved,  but  at  lower 
efficiency  (see  Figure  2.9).  Note  that  these  characteristics  are  all  for 
an  untapered  cavity.  For  a slightly  tapered  cavity,  we  might  expect  an 
efficiency  of  » 30%,  with  an  output  power  of  168  kw. 

Competing  Modes 

One  consideration  not  discussed  above  is  that  of  competing  modes.  All 
of  the  previous  analysis  is  made  on  the  assumption  that  only  the  desired 
mode  will  be  excited.  A mode  of  the  system  exists  only  when  there  is  a near 
intersection  of  the  beam  and  waveguide  dispersion  relations;  i.e.,  when  the 
phase  velocities  are  almost  equal.  For  the  35  GHz  oscillator,  the  disper- 
sion relations  for  the  TEq^  and  neighboring  TEq^^2,3,4. . . “odes  are  given  in 
Figure  2.10.  For  any  of  these  modes  to  oscillate  it  is  necessary  that  the 
threshold  QP^eam  (for  the  mode  of  interest)  to  be  above  that  of  the  system. 
Only  the  analysis  of  azlmuthally  symmetric  modes  has  been  done,  and  thus  the 


Figure  2.8.  Obtainable  and  optimum  Q's  versus  the  ratio  b/L 
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GYROMONOTRON  DESIGN 


Beam  Energy 
Beam  Current 
Frequency 
Cavity  Mode 
Output  Mode 

(without  converter) 
(with  converter) 

Cavity  Q 
Cavity  Length 
Max.  Cavity  Radius 
Efficiency 
Output  Power 

*Value  determined  in  cold  test. 


PARAMETERS 

70  kv 
8 Amps 
35.0  GHz 
TE 

Oil 

“01  (circular) 

TE 

10  (rectangular) 
390* 

2.85  cm 
0.530  cm 
27* 

151  kw 
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Chreshold  condition  can  be  determined  only  for  these  modes.  However,  from 
symmetry  considerations,  It  is  expected  that  the  interaction  for  S-l  will  be 
strongest  for  the  TEQn£  modes.  Additionally,  non-azlmuthally  symmetric  modes 
may  be  preferentially  suppressed  by  interrupting  or  loading  axial  wall 
currents,  which  are  not  present  in  the  TEQn£  modes. 

For  the  35  GHz  oscillator,  the  modes  which  look  most  likely  to  compete 
with  the  TEqh  mode  are  the  TEq^£  modes,  where  1 - 2,3,4...  . The  minimum 
threshold  values  of  QP^,  for  these  modes  are  given  in  Table  2.2.  Also 
given  are  the  threshold  values  for  the  magnetic  field  at  which  the  maximum 
efficiency  is  expected  for  the  TEq^  mode.  The  difference  between  these 
values  results  from  the  fact  that  the  minimum  in  the  threshold  QP^  occurs  at 
a different  phase  shift  for  each  mode.  We  note  that  strong  mode  competition 
would  be  expected  if  the  Q for  all  modes  were  equal.  Fortunately,  Q is 
expected  to  go  as 


where 


and 


f and  f are  the  cutoff  frequencies  of  the  output  guide  and  cavity, 
co  cc 

respectively.  Q is  therefore  a strongly  decreasing  function  of  l,  and  the 
minimum  separation  will  be  Q(J-)  ■ ~ Q(l),  with  p ■ 0.  Therefore,  no  diffi' 
culty  is  expected  from  mode  completion  with  this  oscillator. 


9 


31 


f 


I 


» 


TABLE  2.2 

THRESHOLD  VALUES  FOR  QPb  FOR 
TEQU  MODES 

35  GHz  CAVITY 
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In  leaving  this  subject,  we  note  that  although  the  problem  of  competing 
modes  does  not  seem  serious  for  a gyrotron  operating  In  the  TEq11  mode  with 
s*l,  it  may  prove  to  be  a limiting  factor  in  devices  with  higher  mode 
numbers. 

Mechanical  Design 

The  mechanical  design  of  the  35  GHz  oscillator  was  prepared  with  the  goal 
of  producing  a device  that  was  as  simple  and  flexible  as  possible,  and  that 
could  be  fabricated  in  the  NRL  shops.  The  former  concern  reflects  the 
desire  that  the  cavity  and  other  microwave  structures  (such  as  the  transition 
from  the  cavity  to  the  output  guide)  be  easily  replaceable.  For  this  reason, 
the  microwave  interaction  structures  were  mounted  in  a vacuum  shell,  thus 
removing  the  requirement  that  joints  between  the  structures  be  vacuum  tight. 

The  second  concern  Is  not  trivial.  Vacuums  of  10-8  to  10”  7 torr  are 
required  for  operation  of  the  electron  gun.  The  residual  gases  must  also  be 
free  of  sutstances  which  will  contaminate  the  cathode.  Principal  among  these 
substances  is  zinc,  which  is  a component  of  all  hard  solder  which  can  be  used 
with  a torch.  This  difficulty  can  be  eliminated  by  brazing  materials  in  a 
hydrogen  or  vacuum  oven.  Although  such  facilities  existed  at  NRL,  there 
were  no  trained  personnel  to  operate  them.  Therefore,  it  was  decided  to  make 
all  parts  of  stainless  steel,  and  weld  all  joints.  This  method  has  worked 
reasonably  well  for  laboratory  experiments. 

The  oscillator  is  shown  in  Figure  2.11.  We  shall  discuss  each  component 
in  order,  from  left  to  right. 

Electron  Gun 

The  electron  accelerator  is  of  NEL/Varian  design^ and  manufactured 


by  Varian.  The  design  calls  for  an  operating  voltage  of  70  kv,  with  a beam 
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current  of  8-10A.  Under  these  conditions,  the  ratio  of  perpendicular  to 
parallel  electron  velocities  should  be  1.5,  and  the  spread  in  the  parallel 
and  perpendicular  velocities  to  be  10%  and  7%,  respectively.  Chu^  esti- 
mates that  these  spreads  should  not  significantly  effect  the  maser  mechanism. 

Beam  Input  Guide 

The  guide  diameter  was  chosen  to  be  0.775  cm  (0.305  inches)  in  order  to 
minimize  the  possibility  of  the  ^2^  modes  being  excited.  At  35  GHz,  the 
cutoff  diameter  for  this  mode  is  .802  cm.  In  light  of  experiments  being 
conducted  at  the  time  of  the  writing  of  this  report,  this  may  have  been  an 
unwise  policy,  since  it  appears  that  a substantial  portion  of  the  beam  is 
being  intercepted  by  this  guide.  The  beam  outer  diameter  is  specified  to 

be  .660  cm  (.260  inches),  so  interception  should  not  take  place  with  careful 
alignment.  There  is  a question  then,  as  to  the  actual  size  and  character  of 
the  beam,  and  experlements  are  being  prepared  to  examine  these  parameters. 

Cavity 

Three  cavities  were  fabricated,  two  with  a slight  taper,  one  with  no 
taper.  Their  dimensions  are  given  in  Figure  2.12.  We  note  that  cavity  it 6 
is  the  cavity  specified  in  the  design.  Cavity  it 5 was  chosen  to  test  the  effect 
of  tapering,  and  cavity  it 4 was  chosen  to  test  the  effect  of  varying  L/b. 

The  input  end  of  the  cavities  approximated  a solid  wall,  in  that  the 
input  guide  was  a strongly  cut  off  guide.  The  output  end  of  the  cavity  was 
defined  by  an  outward  step  in  diameter.  This  step  is  a discontinuity  in  the 
guide  impendance,  and  thus  causes  a reflection.  For  the  TE  modes,  the 
impedance  is  lower  in  the  larger  diameter  guide.  Thus  the  discontinuity 
looks  like  a partial  "short  circuit."  The  reflection  caused  by  the  step  is 
given  by  equation  2.10.  This  expression  uses  the  intrisic  impedance  of  a TE 
wave  in  an  infinitely  long  waveguide.  Since  none  of  the  "guides"  are  uniform 
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Figure  2.12.  Diagram  of  the  cavities  used 


in  the  region  of  interest  for  more  than  a few  wave  lengths,  the  expression 
for  the  reflection  coefficient  is  an  approximation  and  must  be  checked  by 
cold  testing. 
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The  Q resulting  from  the  cavity  geometry  is  given  by  2.13.  A Q of ~400 
was  desired.  In  order  to  achieve  this  value  of  Q,  a step  larger  than  that 
provided  by  the  output  waveguide  is  demanded  by  2.13.  However,  cold  tests 
indicated  that  the  step  sizes  given  in  Figure  2.12  be  used.  See  below  for  a 
discussion  of  this  discrepancy. 

Output  Waveguide  and  Collector 

The  output  guide  was  chosen  to  be  a standard  circular  waveguide  size, 
with  a diameter  of  1.603  cm  (0.631  inches).  This  guide  had  approximately 
fifty  0.159  cm  diameter  holes  drilled  in  it  for  pumping.  It  also  serves  as 
a rudimentary  collector  allowing  pulse  rates  of  ~10  pps. 

Output  Window 

Two  configurations  of  output  window  were  designed.  Both  were  one  half  or 
three  halves  wave  length  resonant  windows  fabricated  from  beryllia.  Beryllia 
was  chosen  as  a common  microwave  window  material  with  high  heat  conductivity, 
a property  which  will  be  useful  in  the  long  pulse  gyrotron. 

The  two  windows  are  shown  in  Figure  2.13.  The  first  (Figure  2.13a),  is 
designed  to  be  used  with  the  mode  converter.  With  this  window,  the  output  wave- 
guide is  uniform  in  diameter  from  the  collector  to  the  mode  converter.  The 
second  window  (Figure  2.13b)  is  five  centimeters  in  diameter,  and  was  used  in 
conjunction  with  a horn  for  radiating  the  gyrotron  output.  The  reflectance  of 
both  of  these  windows  was  measured.  The  results  are  shown  in  Figures  2.14 
and  2.15. 
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Figures  2.13a  and  2.13b.  Schematics  of  two  different  output  window 

configurations 
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Figure  2.15.  Plot  of  incident  and  reflected  power  versus  frequency  for 

the  large  window.  The  window  was  mounted  in  a horn  for  the 
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Converter  and  Filter 

The  gyrotron  discussed  in  this  report  was  designed  to  produce  radiation 
in  the  TEq^( circular)  mode.  This  mode,  while  being  very  low  loss,  cannot 
be  focused  and  is  difficult  to  handle.  It  is  therefore  desirable  to  convert 
the  radiation  info  a rectangular  mode,  (TE^q) , which  can  be  focused  and  can 
be  handled  using  readily  available  components.  Two  converters  were  commercially 
available  for  this  purpose.  Their  specifications  are  given  in  Table  2.3. 

Their  reflection  and  transmission  coefficients  were  measured,  and  found  to  be 
within  these  specifications.  The  Hitachi  converter  seems  most  likely  to 
with  stand  high  peak  powers,  and  was  used  in  the  gyrotron  tests.  A mode 
filter,  which  absorbs  all  but  TEQn  (circular)  modes  was  used  to  reduce 
reflections  of  non  TEQn  modes,  if  present. 

TABLE  2.3 

Converter Max  VSWR Max  Conversion  Loss 

Hitachi  1.15  0.5  db 

TRG  1.10  0.3  db 

Detectors-Peak  Power 

Standard  IN53  type  broad  band  crystal  detectors  were  used  for  measure- 
ment of  peak  powers.  These  units  were  individually  calibrated  using  a low 
level  CW  signal  and  thermally  compensated  Hewlett-Packard  thermistor  detector. 

It  should  be  noted  that  the  crystal  detectors  were  used  in  conjunction  with  a 
rectangular  waveguide,  and  thus  measured  only  radiation  originally  in  the  TEq^ 
(circular)  mode. 


Detectors-Average  Power 

A water  calorimeter  was  constructed  for  measurement  of  the  average  gyrotron 
power.  This  device  is  shown  in  Figure  2.16.  Of  particular  Importance  is  the 
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fact  that  the  calorimeter  is  designed  to  absorb  power  in  all  modes.  It  can  thus 
be  used  with  one  of  the  TEQn  filters,  to  determine  (by  subtraction)  the  power 
in  non-azimuthally  symmetric  modes. 

The  calorimeter  consists  of  six  turns  of  plastic  tubing  in  a 1.60  cm 
diameter  waveguide.  Water  was  pumped  through  the  tubing.  A pair  of  thermo- 
couples was  used  to  measure  the  change  in  the  water  temperature  as  it  passed 
through  the  tubing.  The  device  was  calibrated  using  a 20  lew  magnetron, 
pulsed  at  « 100  pps,  and  a themistor  detector.  The  calorimeter  appeared  to 
be  accurate  to  approximately  + 5%.  Care  had  to  be  taken  to  assure  that  the 
water  flow  was  the  same  each  time  the  calorimeter  was  used  in  order  to 
achieve  this  accuracy. 

Test  Results 

Cold  Tests 

Cold  tests  (i.e.,  microwave  measurements  in  the  absence  of  the  electron 
beam)  were  performed  on  the  cavity  and  output  structure. 

The  resonant  frequency,  Q,  and  field  distribution  were  measured  for  the 
two  tapered  cavities.  The  apparatus  used  for  this  measurement  is  shown  in 
Figure  2.17.  Power  was  coupled  into  and  out  of  the  cavity  via  small  (2  mm 
diameter)  holes  on  opposite  sides  of  the  cavity,  near  the  mid-point  of  the 
cavity.  Ka  band  waveguides  were  oriented  to  couple  to  the  H2  of  the  cavity. 

The  resonant  frequency  was  measured  using  a standard  absorbtion  type  wave- 
meter.  The  Q was  determined  by  measuring  3 db  points  of  the  resonance  curve. 
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The  field  distribution  in  the  cavity  was  determined  by  the  perturbation 
method.  In  this  case  a .025  cm  thick  disk  of  plastic,  with  a diameter  almost 
as  large  as  that  of  the  cavity  was  used  as  the  perturbation.  The  local 
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electric  field  magnitude  is  inferred  from  the  change  in  the  resonant  frequency 
due  to  the  presence  of  the  disk.  The  relation  is^7^ 

« . . ±!a  4i 

, £°  2%  / Vdv 

vo 

where  fg  is  the  resonant  frequency  without  the  disk,  Af  is  the  shift  in 

• frequency  from  fQ,  and  £q  are  the  disk  and  vacuum  permitivities , and 

and  Vq  are  the  volumes  of  the  disk  and  cavity  respectively. 

The  cavity  specifications  are  given  in  Table  2.4,  and  the  field  distribu- 

• tion  for  cavity  #4  in  Figure  2.18.  The  resonant  frequencies  are  approximately 
as  expected.  The  Q's  observed  are  not  as  predicted,  as  it  is  expected  that 

V 

the  longer  cavity  should  have  the  higher  Q.  (The  longer  cavity  should  have  a 
Q of  approximately  twice  that  of  the  shorter.)  In  addition,  the  Q's  observed 
are  all  above  that  which  are  predicted  by  2.13.  This  brings  to  question 
either  the  rather  simplified  means  of  calculating  Q,  and/or  the  measurement. 

Both  of  these  possibilities  will  be  pursued,  and  be  reported  on  in  the  reports 
for  project  6085.  For  the  purposes  of  the  initial  experiment,  however,  the 
cold  test  values  will  be  used. 

We  note  that  the  field  profile  for  the  tapered  cavity  is  weighted 
toward  the  output  of  the  cavity;  as  was  expected. 
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Teats  with  the  Electron  Beam  Present 

Preliminary  "hot"  tests  were  conducted  under  this  contract.  Due  to 
the  fact  that  all  of  the  parts  of  the  system,  as  well  as  those  of  the  test 
facility,  were  untried,  most  of  this  testing  was  directed  toward  "debugging." 
However,  the  system  and  test  facility  are  working  well,  and  a full  characteri- 
zation and  optimization  should  be  fairly  straight  forward.  These  results  will 
be  included  in  the  reports  for  JAYCOR  Project  if 6085. 

The  general  behavior  of  the  device  is  as  follows: 

The  device  oscillated  stably,  and  once  the  correct  parameter  regime  was 
determined,  with  very  little  adjustment.  Repeatability  from  run  to  run  was 
excellent,  demonstrating  the  better  than  one  percent  stabilities  of  the 
power  supply,  modulator,  and  magnet  controller.  Under  optimized  conditions, 
the  output  pulse  was  quite  flat  and,  within  the  admittedly  poor  resolution 
(~  10  MHz)  of  the  cavity  wavemeter,  fairly  monochromatic.  A better  measure- 
ment of  the  spectrum  was  attempted  using  a spectrum  analyzer,  but  the  results 
obtained  were  Inconclusive,  probably  due  to  Instrumentation  difficulties. 

The  two  cavities  that  are  shown  in  Figure  2.12  were  tried.  The  data 
obtained  with  each  was  not  complete  enought  to  determine  if  the  difference 
in  lengths  is  significant.  A cavity  with  no  taper,  but  with  a length  the 
same  as  the  shorter  tapered  cavity  was  fabricated,  but  not  tried. 

The  maximum  power  obtained  was  87  kw;  with  an  efficiency  of  12%.  The 
maximum  efficiency  obtained  was  13Z,  at  77  kw.  The  pulse  length  was  approxi- 
mately 1.5  /usee,  and  the  repetition  rate  £ 10  pps.  The  frequency  of 
oscillation  was  35.05  GHz,  with  a tunability  of  approximately  + 150  MHz. 

The  efficiencies  and  output  powers  are  clearly  disappointing.  However, 
measurements  made  on  a long  pulse  version  of  the  oscillator  Indicate  that 


up  Co  40  percent  of  Che  electron  beam  Is  being  intercepted  at  the  Input  guide. 
This  will  clearly  alter  the  true  efficiency,  and  raises  serious  questions  as 
to  the  actual  beam  size,  position,  and  quality,  since  by  design  the  beam  should 
have  little  difficulty  passing  through  the  input  guide.  At  the  writing  of  this 
report,  tests  are  in  progress  to  determine  the  exact  nature  of  the  problem. 

The  results  of  these  tests,  and  those  of  measurements  made  with  knowledge  of 
the  beam  current  will  be  detailed  in  subsequent  reports. 
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Conclusion 

Under  contract  No.  N00173-77-C-0215,  JAYCOR  has  developed  methods  for 
the  design  of  gyrooonotrons  with  right  circular  cavities.  Using  these  methods, 
a gyromonotron  has  been  designed  to  operate  at  35  GHz  with  an  output  power 
of  170  kw  with  an  efficiency  of  302.  This  efficiency  is  in  rough  agreement 
with  those  predicted  by  others^*^^,  for  similar  devices. 


1 


The  gyrotron  has  been  fabricated  and  a gyrotron  test  facility  designed 
and  constructed.  The  device  has  undergone  initial  tests.  The  results  of 
these  tests  have  indicated  much  lower  efficiencies  than  predicted,  but  are 
most  likely  complicated  by  the  interception  of  the  beam  prior  to  entry  into 
the  cavity.  It  is  believed  that  this  problem  can  be  corrected  fairly  easily, 
and  complete  results  of  subsequent  tests  will  be  reported  under  JAYCOR 

Project  #6085. 
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APPENDIX  A 

MASER  TESTING  FACILITY 

A maser  test  facility  has  been  established  and  tested  at  the  NRL  site. 

The  basic  system  consists  of:  (a)  a 5 Mw  lumped-line  modulator  and  pulse 
transformer  combination,  which  has  been  loaded  and  tapped  to  give  70  kv@ 

0-15A  and  35  kv  + 3500  v @ 1A  with  a fixed  pulse  length  of  approximately 

2 ps,  (b)  a 40  kg  superconducting  magnet.  The  magnet  specifications  are  on 

file  in  room  104  of  building  51  at  NRL,  and  the  approximate  field  profile 

is  given  in  Figure  2.11,  (c)  auxiliary  subsystems,  including  timing  circuits, 

a heater  current  source  and  diagnostics,  etc.  These  items  are  conventional, 

and  their  circuit  schematics  are  on  file  in  room  104  of  building  51  of  the 

Naval  Research  Laboratory.  y 

The  facility  is  presently  fully  operational.  It  is  a straightforward 
system  which  has  required  little  maintenance,  and  should  provide  good  continued 
service. 

Manuals  for  the  individual  components  of  the  sytem  are  located  with 
those  instruments,  in  room  104  of  building  51  at  the  Naval  Research  Laboratory. 


51 


REFERENCES 


1.  P.  Sprangle  and  A.  T.  Drobot.  "The  Linear  and  Self-Consistent  Non- 
linear Theory  of  the  Electron  Cyclotron  Maser  Instability";  IEE  Trans. 
MTT-25 . 528-544  (1977). 


2.  K.  R.  Chu.  "Theory  of  Electron  Cyclotron  Maser  Interactions  in  a 

Cavity  at  the  Harmonic  Frequencies",  NRL  Memo  Report  3672,  (1977). 


3.  V.  A.  Flyagin,  A.  V.  Gaponov,  M.  I.  Pelelin,  and  V.  K.  Yulpatov. 

"The  Gyrotron",  IEE  WTT,  25,  1977. 

4.  S.  A.  Schelkunoff.  Electromagnetic  Waves,  Van  Nostrand  Co  (NY),  1943. 

5.  J.  L.  Seftor,  K.  R.  Chu,  and  A.  T.  Drobot.  "An  Investigation  of  a 

Magnetron  Injection  Gun  Suitable  for  Use  in  Cyclotron  Resonance  Masers", 
NRL  Memo  Report  3697,  (1978). 

6.  Varian  Associate,  "Technical  Report  on  the  Design  of  an  Electron  Gun  for 

a Cyclotron  Resonance  Maser",  under  NRL  Contract  N00173-77-C-0086. 

(June  1967). 


7.  M.  A.  Heald  and  C.  B.  Wharton,  Plasma  Diagnostics  with  Microwaves. 
John  Wiley  and  Sons,  New  York,  1965. 


